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Abstract 



Recently it was pointed out that the solutions found in literature for the 

space fractional Schrodinger equation in a piecewise manner are wrong, except 

the case with the delta potential. We reanalyze this problem and show that an 

exact and a proper treatment of the relevant integral proves otherwise. We also 

K^ i discuss effective potential approach and present a free particle solution for the 

>0 ' space and time fractional Schrodinger equation in general coordinates in terms 

1^ . of Fox's H-functions. 
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I. Introduction 

In a recent article, Jeng et. al. [1] argued that solutions obtained for the 
space fractional Schrodinger equation in a piecewise fashion are not valid. Their 
arguments were based on a contradiction that they observed when the ground 
state wave function of the infinite potential well problem is substituted back into 
the space fractional Schrodinger equation. They concluded that many exact so- 
lutions presented in literature are wrong, except the one with the delta function 
potential. This not only casts doubt on the existing solutions in literature, but 
also makes it very difhcult to find a meaningful solution. In this manuscript, 
we show that an exact treatment of the integral that lead them to inconsistency 
proves otherwise. 

Fractional calculus has found to be extremely useful in the study of anoma- 
lous diffusion. The fastest way to introduce memory effects into an existing 
equation is to replace a time derivative with its fractional counterpart. In this 
regard, recently we have reconsidered the time fractional Schrodinger equation 
and its separable solutions [2]. The Space fractional Schrodinger equation nat- 
urally incorporates non local effects and was first introduced by Laskin [3-7]. 
It is intriguing that the path integral formulation of quantum mechanics over 
Levy paths, which has many interesting examples in nature, leads naturally to 
the space fractional Schrodinger equation where the space derivative is replaced 
by the fractional Riesz derivative [2,7]. In this regard, Riesz derivative plays a 
central role in the fractional generalization of the Schrodinger equation, where 
solutions for various potentials are available [3-9] . It should be emphasized that 
even thought memory and/or non local effects can be studied via the fractional 
generalizations of certain derivatives in physical equations, there could also be 
interesting cases where fractional derivatives alone are not sufficient [10,11]. 
For example, path integral formulations of quantum mechanics over distribu- 
tions other than the Levy distribution, would not necessarily lead to a fractional 
Schrodinger equation. Fractional extension of quantum mechanics that treats 
coordinates and conjugated momentum equivalently have also been considered 
in literature [12]. 

To understand the effects of fractional derivatives in Schrodinger equation, 
one could try effective potentials within the framework of quantum mechan- 
ics. In this regard, we have introduced the effective potential approach for the 
time fractional Schrodinger equation [2], which allowed us to write the time 
dependence of the separable solutions, which is a Mittag-Leffier function with 
an imaginary argument, as the product of its purely oscillating and purely de- 
caying parts. In this manuscript, we also write the effective potential for the 
space fractional Schrodinger equation and apply it to the infinite well problem. 
We also present a general solution for the space and time fractional Schrodinger 
equation in terms of Fox's H-functions in general coordinates. 



II. Consistency of the Solutions of the Space Fractional 
Schrodinger Equation 

It is well known that the one dimensional space fractional Schrodinger equa- 
tion is given as 

^n^^D„{-h'Ay^'^ + V^, (1) 

where (—h^A) is the quantum Riesz derivative [7,13]: 

1 r-\-oo 

(-fi2A)"/'vl/(a;,t) = _y dpe^P-/>'\pr^p,t), (2) 

and ^{p,t) is the Fourier transform of '^{x,t) : 

/ + 00 
dx^{x,t)e''P'=/''. (3) 

-oo 

For the infinite square well potential: 

r if |x| < a 

V{x) ^1 , (4) 

[ oo if \x\ ^ a 

the separable solutions are of the form 

*(x,t) = e-*^*/V(^), (5) 

where tp{x) satisfies the eigenvalue problem 

Da {-h^A)"^^ tl^ix) = Eil^ix), V(a) = V(-a) = 0. (6) 

The general solution is obtained as [4,7,8] 

Acos{^{x + a)) if |x|<a 

M^) = { " , n = l,2,... (7) 

if \x\ ^ a 



with the energy eigenvalues 



©°- <«' 



To show the inconsistency of these solutions, Jeng et. al. [1] concentrated on 
the ground state with n = 1: 

ylcos(-l^) if \x\ < a 

Mx) ={ . (9) 

if \x\ ^ a 



and argued that this solution, albeit satisfying the boundary conditions, V'l (—a) = 
tpi{a) = 0, when substituted back into the space fractional Schrodinger equation 
[Eq. (6)] leads to a contradiction. 

The Fourier transform of ipi (x) can be found as 



fi^ \ cos (ap/H) 
a ) p^- {TTn/2af 



dp)--A.(-] . :"r:\, , \x\<a, (10) 



which when substituted into Equation (6), gives ^'1(2^) as the integral 

Defining a new variable, q ~ -^P, we can also write this as 

V,,(:,) . -^£^ f^Y r°^dq M:£^W^^..W2a_ (12) 

TiEi \2aJ J_^ q^ -1 

Without evaluating the above integral, Jeng et. al. [1] argued that the right 
hand side, which they wrote as 

ttEi \2aJ Jo '^ q^-l 



M-) = -^i^] 2/ dg '^' 7^'7-^ cos(.W2«), (13) 



can not satisfy the boundary conditions as a; — > ±a that "01(2;) satisfies, thus 
pointing to an inconsistency. However, an exact evaluation of the integral in 
Equation (12) proves otherwise. Note that the integral 

J -00 q - 1 

is a singular integral with poles on the real axis at q — ±1, which could be 
evaluated via analytic continuation as a Cauchy principal value integral [13]. 
Substituting 

cos (TTq/2) = i U'"'/^ + e-'"""/^) , (15) 

we first write / as two integrals: 

= /i+/2. (17) 

For the first integral, /i, we close the contour in the upper half q— plane with a 
semicircular path of radius R and then go around the poles on the real axis in 
the upper q— plane, with small semicircular paths of radius S. In the limit as 
i? — >■ 00, via Jordan's lemma, the contribution over the large semicircle vanishes. 



thus allowing us to evaluate the value of this integral as a Cauchy principal value 
in the limit (5 ^ as [13] 

Pm)--Qcos|^). (18) 

Similarly for I2, we close the contour this time in the lower g— plane and circle 
around the poles in the lower half g— plane to obtain 

m/.)--(4cos^). (19) 

Adding these we obtain the Cauchy principal value of / as 

TTT 

PT/(/)--^cos— , (20) 

Za 

which on the contrary to Jeng et. al. [1] vanishes at the boundary as x — > ±a. 
When this is substituted back into Equation (12) yields 

^i(x) = ^-(^-j cos-, N<a. (21) 

Since the eigenvalue for the ground state is 

E.-D.[^-) , (22) 

we again obtain the ground state wave function [Eq. (9)]: 

Tpi{x) ^ Acos—, |a;| < a, (23) 

which satisfies the boundary conditions, hence no inconsistency. 

The proof for general n also follows along the same lines. For example, for 
the odd values of n. Equation (12) is replaced by 



Now the Cauchy principal value of the needed integral: 

'•+°° \q\"cos{mrq/2) j(^^^/2a)c, 

is found as 



/ = / dq 1^1 — V-^'^/-y ^.(nW2ak^ (25) 



m/) = — m(-)cos(— ), (26) 

which when substituted back into Equation (24) yields 

, , , ADa . o /n7r\ fmrH\ /mrx\ . . 

^-j, n=l,3,5,.... (28) 

Again, no contradiction. The proof for the even values of n is similar. 



III. Effective Potential for the Space Fractional Schrodinger 
Equation 

For the separable solutions of the free space fractional Schrodinger equation 
we use Laskin's notation: 

ift^*(a;, t) - -i?^(W)%(x, t), (29) 

where V^ is the Riesz derivative and 

(-fi^A)"^' ^{x, t) = -(ftV)'3*(x, t). (30) 

For a separable solution, ^(x,i) = T{t)X{x), found by solving the free space 
fractional Schrodinger equation we can define an effective potential in terms of 
the Schrodinger equation as 

h^ I ( dPX{x) \ D,{hVYX{x) 
^^"■^""^ ^ 2^X{x) [^^) X{^) ■ (^') 

For separable solutions, the free space fractional Schrodinger equation reduces 
to 



ih^^ - Er,T{t), (32) 

-DfiihVfXix) ^ EnX{x), (33) 

where the energy eigenvalues, _E„, come from the solution of Equation (33) with 
the appropriate boundary conditions and, the time dependence is given as 

T{t) = e-'(-^"/'^)*. (34) 

Since in the presence of a potential, V{x), the corresponding fractional Haniil- 
tonian: 

H^^^DfiihV)^ + V{x), (35) 

is hermitian, the energy eigenvalues are real [7]. Thus, we can also write 
The equivalent Schrodinger equations are now written as 



dt 
1 (dPX{x 



= Er,T{t), (37) 



2m X{x) \ dx"^ 
which are to be solved with the same boundary conditions 
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+ Veff.{x)^Er„ (38) 



A. Effective Potential for the Infinite Well 

For the box problem, using the normaUzation condition 



J — a 


= 1, 




(39) 


the complete wave function is given as 








r 1 e-(V^)^"*sin7(x + a) 


; 


\x\ < a 


,n = l,2, 





1 


\x\ > a 


(40) 



where 



mrn 



^" " ^^ ' ^7 ' ''^'^ *^""' ^^ " *^"' *^ " °- ^^^^ 



The effective potential is now written as 

In other words, the effective potential shifts the energy levels by a constant as 
a function of the quantum number n. 



IV. Free Particle Solution of the Space and Time Frac- 
tional Schrodinger Equation 

Wc now consider the free Schrodinger equation with both time and space 
fractional derivatives [2]: 

^*(x, t) = '-j-Dcji {hVf ^(x, i), < a < 1, 1< /3 < 2, (43) 

where 4^ is the Caputo derivative (Appendix A). Performing a Wick rotation 
gives the Bloch equation: 

^*(x, t) = ^b^.ph^S/i^ix, t), < a < 1, 1< /3 < 2, (44) 

where D\^2 = 1/2to. Using the boundary conditions 

* (x, 0) = (5(a;) and lim *(x,t)^0, (45) 

a;— >±oo 

we take the Laplace transform with respect to time and the Fourier transform 
with respect to space to obtain the Fourier-Laplace transform of the solution 



[14]. Finding the inverse transform and then performing an inverse Wick rota- 
tion, yields the wave function in integral form as 



*(x,i) 



+ 00 



J r''e~''''Ej-'-D^.,hPk^e] dk, 



(46) 



where Ea{z) is the Mittag-Lcffler function [2]. We can also write 'ii{x,t) as 
*(x,i) == — / cos /ex Ea, -—Da,ph^kH'^ dk. 



(47) 



This wave function satisfies both the time and the space fractional Schrodingcr 
equation [Eq. (43)] 

In terms of iJ— functions (Appendix B) Equation (47) can be written as 



+00 



^{x,t) ^^ I coskx hI'^ ( "--Do^^ph^kh" 



(0,1) 



dk. 



(48) 



"'O \ " (0,1), (0, a)/ 

which can be integrated by using the properties of the H— functions [15] as 



M..,:)^^^Rli[';b...n'eil 



(l/2,/3/2),(0,l),(0,/3/2) 



(0,1), (0,a) 

To determine "So the final wave function has to be normalized as 



/|^(x,0)|^.x^l. 



(49) 



(50) 



Case I: 

For /3 = 2, this solution becomes the free particle solution of the time frac- 
tional Schrodingcr equation as 



^ r ^ *0 ,a 2 / 4i"i:)„t" 

*(x,t) = . H ' 



'3,2 



(1/2,1), (0,1), (0,1)N 
(0,1), (0,a) 



Da = Da^fl 



(51) 



This can also be shown to be equal to 
^(x,i) 



*0 jj2,0 

7^TZT^i.2 



Ai'^Dat" 



or to 






(1,") 

(1/2,1), (l,l)y 
(l,a)N 

(1,2), 



(52) 



(53) 



In the limit as a — > 1, this becomes 



*'•'•'> ^iiJio^'^^-'f-in^)- f^"' 



where Di — h/2m. 
Case II: 
When a = 1, we obtain the wave function as 



*n „i W «„.«./ 2 



*<-')^iH''a >nR 



/3 



(l/2,/3/2).(0,l):(0,/3/2)N 



(0,1), (0,1) 



(55) 



which satisfies 



where 



i h—^{x,t) = -DfjH>^R^^{x,t), 1 < /3 < 2, (56) 



Di,p = Df,. (57) 

This solution was also given by [8] in the form: 

(1,1//3),(1,1/2)N 



TT^o ^14 / 1 / h ^ '/^ 



p\x\ ' \ h \iDfit 



(58) 



(1,1), (1,1/2) 



V. Conclusions 



Successful applications of fractional calculus to anomalous diffusion attracted 
researchers from many different branches of science and engineering into this 
intriguing branch of mathematics [15-19]. Applications of fractional calculus 
usually starts by replacing certain derivatives in the evolution or the transport 
equations with their fractional counterparts. In general, replacing a time deriva- 
tive with its fractional counterpart incorporates memory effects into the system 
and makes the process non Markovian, while a replacement of a space derivative 
introduces global or non local effects. 

In 2000, Laskin [3-7] introduced the path integral formulation of quantum 
mechanics over Levy paths and showed that the corresponding equation of mo- 
tion is the space fractional Schrodinger equation. For the infinite potential well 
problem, the complete wave function is given as [7,8] 

-= sin — (x + a) , jxj < a 

*(x,t)=<( Va 2a ,n = l,2,, 

, |x| > a 

(59) 



Since the energy operator is now given as ih—-, the energy eigenvalues are given 

at 
as 

^"^^^(?) ' " = 1'2'---- (60) 

Note that for the space fractional Schrodinger equation, the Hamiltonian, H = 
—Dfi (ft^) + V{x), is hermitian. An alternate approach to space fractional 
Schrodinger equation is given by Herrmann [12]. 

It is true that the Riesz derivative requires a knowledge of the wave function 
over the entire space. For the infinite well problem, the system is confined to 
the region |a;| < a with Vl/(a;,t) = for \x\ > a. Since the solution for |x| < a, 
also satisfies the boundary conditions as a; — > ±a, the solution inside the well is 
consistent with the outside. In this regard, it is not true that the above solution 
ignores non locality and is valid only for /? « 1 [1,12 pg. 107]. 

In general, in the presence of a potential, V{x), the time independent space 
fractional Schrodinger equation is given as 

Dc,{-n^A)°'^^i:{x;E,a) + V{x)i:{x; E,a) ^E^{x]E,a), (61) 

which is to be solved with the appropriate boundary conditions. The Fourier 
transform of Equation (61) yields a dispersion relation of the form (f>(j);E,a), 
where the wave function is given as 

1 r^°° 

Mx; E,a) = ^-^J dp 0(p; E, a)e'P-/\ (62) 

Since quantization follows from the boundary conditions that ')p{x;E,a) satis- 
fies, granted that the above integral exists, solutions satisfying the given bound- 
ary conditions can not lead to an inconsistency [7-9] . 

After showing the consistency of the solutions of the space fractional Schrodinger 
equation for the one dimensional infinite square well, we have also given the 
effective potential for the space fractional Schrodinger equation. This is the 
potential that the Schrodinger equation with the same boundary conditions will 
yield the same wave function as the fractional case. For the space and time 
fractional Schrodinger equation, we give the free particle solution in terms of 
//—functions, which in the appropriate limits reproduces the previous solutions. 
The solutions we give for the free particle are in general coordinates. 
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APPENDIX 

Appendix A. Basic Definitions of the Fractional Deriva- 
tives and Integrals 

The Caputo definition of fractional derivative is given as 



dti 



1 /■* fdf(T)\ dr ,, , 

I I ■'^ ' ^ 0<g< 1, (Al) 



r(i -q)Jo V dT ) [t- t) 



Laplace transform of the Caputo derivative is 

fe=0 



, n-1 < q <n, (A2) 

t=o 



where outfit) = -^• 

In writing the space fractional diffusion equation or the space fractional 
Schrodinger equation, we use the Riesz derivative which is defined with respect 
to its Fourier transform: 

J-{R^/(x)} = - H^icj), 0<q<2. (A3) 

This yields the Riesz derivative as 

n,lf{x)^-—J^ \co\''g{oj)e^-^doj, (A4) 

where 5(0;) is the Fourier transform of f{x). Note that 

Rl/(^) - ^fi^)- (A5) 

Appendix B. Fox's H-Function 

In 1961 Fox introduced the _ff— function, which gives a general way of ex- 
pressing a wide class of functions encountered in applied mathematics. iJ— function 
provides an elegant and an efficient formalism to handle problems in fractional 
calculus. Fox's _ff— function is a generalization of the Meijer's G— function and 
is defined with respect to a Mellin-Barnes type integral [15]: 

Trm.n/ \_ TTm,n f \(.ap,Ap)\ _ „rn,n f \{ai,Ai),...,{ap,Ap)\ (JiT\ 

^p,q y^) - J^p,q \^\(bg,B,)J ~ -"p,q \^\(bi,Bi),...,(b„Bq) J \^^) 

= ^J h{s)z-'ds, (B2) 



11 



where 



his) 



n;ii r(5, + B,s) uu r(i - % - ^.-5) 



m 



B,s) 



(B3) 



m,n,p,q are positive integers satisfying 0<n<p, l<TO<q, and empty 
products are taken as unity. Also, Aj, j = 1, . . . ,p, and Bj, j — 1, . . . , g, are 



real positive numbers, and aj, j 
complex numbers satisfying 



l,...,_p, and bj, j = 1, 



, q, are in general 



^j{i'h + ^) 7^ Bfi(aj — A — 1) for i^, A = 0, 1, . . . ; h = 1, . . . ,m, j — 1, . . . ,n. 

(B4) 
The contour C is such that the poles of T{bj + Bjs), j ~ I, . . . ,m, are separated 
from the poles oir(l — aj — Ajs), j = 1, . . . ,n such that the poles of r(5j +Bjs) 
lie to the left of C, while the poles of r(l — gj — Ajs) are to the right of C. The 
poles of the integrand are assumed to be simple. The _ff— function is an analytic 
function of z for every \z\ ^ when /i > and for < |z| < 1//? when /i = 0, 
where [i and /? are defined as 



M 



E^^- 



p 



(B5) 



and 

Fox's _ff— function is very useful in the study of stochastic processes and 
in solving fractional diffusion equations (Glockle and Nonnenmacher [20] , West 
and Grigolini [21]). For example, the following useful formula for the Riemann- 
Liouville fractional derivative of the iJ— function: 



R-LY>,li 







D 



a Tjm.n 



z^H, 



(cz) 



fj,(aj,Aj) 
kbj.B,) 



^a-P TTm,n+l 



= Z 



h: 



P+1,9+1 



[cz 



u,{-a,b),(aj,Aj) 



\(b, 



,).{P-a,b) 



(B7) 

where a, 6 > and a + bm.va[bj/Bj) > —1, 1 < j < m, can be used to find 
solutions to the fractional diffusion equation by tuning the indices to appropriate 
values. Similarly, the Laplace transform of the if —function can be obtained by 
using the formula [15] 



■^ \^ "p.q+1 {"-^ I(b,,i3,),(l-p 


W _ -p-firn.n ( „„-<y i°-P'^p)\ 
a))]-'' ^P,q l^O'* (6„B,)J' 


(B8) 


where the inverse transform is given as 




-f- \S ^p,q {^^ \(b,,B,) J J - ^ ^p+l,q[aX ^^^^^^^ j, 


(B9) 


where 




p,a,s e C, Re(s) > 0, ct > 


(BIO) 


and 




Re(/9) + (J max 

l<i<n 


■ 1 Re(a,)- 
Ai M 


> 0, arga < — , 9 — a — a. 


(Bll) 
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